Effect of bulk inversion asymmetry on the Datta-Das transistor 
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A model of the Datta-Das spin field-effect transistor is presented which, in addition to the Rashba 
interaction, takes into account the influence of bulk inversion asymmetry of zinc-blende semiconduc- 
tors. In the presence of bulk inversion asymmetry, the conductance is found to depend significantly 
on the crystallographic orientation of the channel. We determine the channel direction optimal for 
the observation of the Datta-Das effect in GaAs and InAs-based devices. 
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Various spin-orbit (SO) effects which are always 
present in semiconductor quantum structures, are recog- 
nized as an important issue in the emerging field of spin- 
tronic devices. Already the first proposal of a spin tran- 
sistor by Datta and Dasi relied on the SO-induced split- 
ting of the carrier bands in high electron mobility transis- 
tors (HEMT). As postulated by Bychkov and Rashba^ 
in such systems, owing to the SO interaction and interfa- 
cial electric field E (the structure inversion asymmetry - 
SIA), the spins of the electrons fiowing with the velocity 
V precess in a magnetic field Bji ^ E x v. To take ad- 
vantage of this fact, Datta and Das considered a hybrid 
structure consisting of a ballistic quasi-one-dimensional 
(ID) wire patterned of a 2D HEMT and of ferromagnetic 
source and drain ohmic contacts, serving as spin injector 
and detector, respectively. The contacts were assumed 
to be magnetized along the current direction and to have 
100% spin injection/detection efficiency. The calculated 
conductance for such a device was found to oscillate as 
a function of the SO coupling strength a, whose magni- 
tude could be controlled via a metal gate located on the 
top of the structure. Datta and Das showed that in the 
lowest order in a the oscillations depended on neither ID 
sub-band index n nor electron wave vector k. 

Since the original proposal of Datta and Das, the oper- 
ation of their spin field-effect transistor (spin-FET) has 
been analyzed in much detail. For example, calculations 
of device conductance not only for small but also for large 
a were carried out^ as well as possible effects of electron- 
electron interactions were analyzedi^ Furthermore, an 
important problem of spin injection across the ferromag- 
net / semiconductor interface was considered from various 
viewpoints However, in addition to the SIA (Rashba) 
term, owing to the lack of inversion symmetry in bulk ma- 
terials - there exists the so-called bulk inversion asym- 
metry (BIA) or Dresselhaus termrS It is known that the 
combined effect of both terms modifies considerably spin 
properties of zinc-blende heterostructuresi^ Therefore, 



the SIA and BIA contributions should be considered on 
equal footing when analyzing the performance of spin 
devicesfi^ 

In this work, we consider a ballistic quantum wire pat- 
terned of a zinc-blende heterostructure grown along the 
[001] direction. The Hamiltonian describing the electron 
motion takes into account both SIA and BIA terms. By a 
numerical solution of the Schrodinger equation for GaAs- 
and InAs-based systems, we demonstrate the importance, 
of hitherto neglected, non-linear terms in the BIA hamil- 
tonian. Moreover, we show that the device performance 
depends strongly on the wire direction in respect to the 
crystal axes. This allows us to determine channel orien- 
tation corresponding to the largest transconductance for 
both GaAs and InAs based devices. 

In order to model the Datta-Das transistor, we con- 
sider an infinite quantum wire that extends in the x di- 
rection, which is oriented at an angle of (p to the [100] 
crystal axis. The potential confining electrons in y direc- 
tion, V{y), is assumed to be parabolic and characterized 
by the energy difference between consecutive harmonic 
oscillator levels hLJo=l meV in GaAs and fiLUo=3 meV 
in InAs, respectively. In order to simulate the presence 
of ferromagnetic electrodes, we assume that in the re- 
gions a; < and x > L the electron spins experience the 
magnetic field |i?out| = lO'^ T, which is strong enough to 
ensure the entire spin polarization. Following Datta and 
Das we take the vectors Bo„t to point in the x direction 
and consider the case of their parallel or anti-parallel rel- 
ative orientations. In the channel region, < a; < L, 
Bout decays adiabatically to zero. For x \\ [100], the SO 
Hamiltonian is given hy^ 



Hso = ia {(Jydx - a^dy) 

-^l{<^A{dl + {kl))-aydy{{kl)-dl)), 



(1) 



where (Jx and ay are the PauH matrices. The terms pro- 
portional to a and 7 corresponds to the SIA and BIA, 
respectively, whose effect on the band splitting has been 
thoroughly examinedii^ Here, we estimate the magnitude 
of a and (fc^) according to Ref . ITll 
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{kl) = i(16.5^e2miV,/Kn2)i/3^ (3) 

which gives values comparable with those obtained from 
the more refined theories, at least in low Fermi energy 
limit. In Eq. k is the static dielectric constant, m is 
the efi'ective mass, Eg is the energy gap, and A is the 
spin-orbit splitting in the valence bandii^ For example, 
for typical areal electron density iVg = 5 x 10^^ cm~^, this 
approach implies a —1.7 meVA and 42.0 meVA for GaAs- 
and In As-based heterostructures, respectively. The BIA 
parameter 7 is equal to 25 eVA'^ for GaAs and 130 eVA3 
for InAspiLii 

By transforming the Hamiltonian (QJ to a new coordi- 
nate system such that the angle between new and old x 
axes is (h we obtain, 



Hso = {aydx - cr^dy) 



-«7 



<Ty sin(20) -a 



1-12 



{kl) 



Gx cos 



{2^) {dl + {kl)) 



-17 



Oy cos 



{dl + {kl)) 



Ox sin(20) ( -d'y 



:dl + {kl) 



a, 
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i7sin(2(?!))fc^(n)a:rX« where kl{n) = 2™ [Ep - e„) is the 
wave vector for n'th mode in the absence of spin-orbit 
coupHng. This approximation is justified by the fact that 

U3 ^"4 Jij. 



for typical experimental situations, ^kp « 
deed, for GaAs for kp ~ 10*m~^ and 7 = 25 meVA^, 

^1^3 _ . A n ^ n-24 



h-'ki 



p ~ 4.0 X 10 J, which is much smaller than „ 
X 10^^^ J. In other words, the spatial changes of the 



wave function are governed mainly by the 



in the Hamiltonian. Treating x„ 
vector X we arrive to 



term 

as the components of a 



-—R{dlx + 2iQdxx) + {h-Ep)x = 0, (8) 
where matrix Q = R^^P and 
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<7y cos(2(/)) — —Ox sin(20) 



Pun' = -^i{5n7i'(JySm{2(j})[{kl) - -kl{n)] 
~S„,i'{kl)<7x cos(20) - [(Tx cos{2<j)) 



-aySm{2^)]Dl^,} 



nn' ' 



(9) 
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We see that only the Dresselhaus (BIA) term changes, 
while the Rashba (SIA) term remains invariant with re- 
spect to the rotations around the z axis. 
The full Schrodinger equation for our calculations reads. 



-^idl + d^+V{y)+Hso 



g^lBCT ■ B„^,{x) *(a;,y) ^ Ep^!{x,y), 



(5) 



where Ep is the Fermi energy. 

To calculate conductance in the non-interacting elec- 
tron model we employ the recursive Green function 
methodic in a form developed in Ref. 0. We expand 
the spinor wave function ^{x,y) into a series 



'^{x,y) ^^X7i{x)'4)7i{y), 



(6) 



where ip7i{y) are the eigenfunctions of one dimensional 
harmonic oscillator potential 



(7) 
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By using the orthonormality property of these functions, 
it is easy to derive a set of differential equations for X7i{x)- 
The problem we encounter is the appearance of 
third order derivative with respect to x in the term 
— 17 sin(20)a^X„(a;), which originates from the BIA 
part of the Hamiltonian. We replace this term by 



hnn' = ( + ■^gfJ'BO- ■ S.,„,(a;) ) 5„ 



—17 



{fc,2)(a,cos(20)-f fT.sin(20))D5 



-^axSm{2q^)Dl^l 



and the the matrix Dmn = {x77i\9y^'\Xn)- After dis- 
cretization along X direction (with the lattice constant 
a) we finally obtain 

- tVx,x+aX{x + a) - tUx,x-aX{x - a) + 2tRx{x) 



iacJxD 



(1) 



(11) 



2m 



RQ\{x) + {h-EF)x{x)^Q, (12) 



where t = Vx,x+a = i?e'"'?, Ux,x-a = Re-'^^Q . 

This form is appropriate for the conductance calculation. 

Before presenting our results we discuss briefly the 
"linear" approximation to the BIA term. Within this 
approximation second order derivatives and dy are 
neglected in Eq. Q . On the basis of this Hnear approx- 
imation one can easily understand the influence of the 
channel crystallographic orientation on the performance 
of spin-FET. Let us consider a transversal running mode 
characterized by wave vector k. Following Datta and 
Dasi we neglect intermode coupling caused by spin-orbit 
terms. Then, it is possible to deflne Rashba and Dres- 
selhaus "magnetic fields" = Q;fc(0, — 1) and Bp, = 
7(A:^)fc(-cos20,sin20). Their sum Bso = Bji + Bp, 
determines the spin precession frequency. Importantly, 
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FIG. 1: Source-drain conductance of a GaAs-based device as 
a function of the Fermi energy for <^=0 ([100] channel direc- 
tion) and anti-parallel electrode magnetizations. Upper panel 
- comparison of linear and nonlinear cases. Lower panel - par- 
tial contributions of three chosen modes to the conductance. 
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FIG. 3: Source-drain conductance of a InAs-based device as a 
function of the Fermi energy for = 0° and = 135° for the 
opposite (upper panel) and parallel (lower panel) source-drain 
polarizations. 
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FIG. 2: Source-drain conductance of a GaAs-based device as 
a function of the Fermi energy for (p = 135° (upper panel) 
and for for (j> = 110° (lower panel). 

the magnitude and orientation of Bso with respect to 
the X axis depend on the angle (j). That is why the pre- 
cession frequency varies with the channel crystallographic 
orientation. It is easy to find out that for (j) = 135°, Bso 
is perpendicular to the x axis, and its length assumes a 
maximum value. Thus, for the [110] channel orientation, 
in the Hnear approximation, the period of the oscillations 
is the shortest, and thus the Datta-Das effect best visible. 
These qualitative considerations have been confirmed by 
our numerical calculations performed in the linear ap- 
proximation. However, in neither GaAs nor InAs the 
linear approximation is justified for the experimentally 
relevant Fermi energy range. 

We start the discussion of the full nonlinear model from 
results for the device oriented along the [100] direction 
{(/) = 0). As shown in the upper panel of Fig. 1, for a 
given relative magnetization of the contacts, two types 
of conductance oscillations as a function of the Fermi en- 
ergy appear: the long-period Datta-Das oscillations and 



the short-period Fabry-Perot oscillations resulting from 
interference between transmitted and reflected spin wave 
functions by the two contactsii In real experiments, the 
latter are expected to be washed out by thermal broad- 
ening and effects of non perfect geometry. Therefore, in 
the following pictures we show smoothed data obtained 
by averaging with the Fermi energy window of 1.0 meV. 
We realize that in the nonlinear case, contributions to 
the conductance originating from particular modes oscil- 
late with different periods. When the Fermi energy or, 
equivalently, electron concentration grows, every 1 meV 
a new mode, with its own period, starts to participate 
in the transport. We illustrate this effect in the lower 
panel of Fig. where a partial conductance due to 0th, 
6th and 12th mode is depicted. This leads to the disap- 
pearance of the conductance minima when the nonlinear 
terms in the SO Hamiltonian are taken into account, as 
shown in Fig. |2l (upper panel). By analyzing the con- 
ductance for different angles (f) we have found that the 
optimal angle for the observation of Datta-Das effect in 
GaAs is (/) = 110° (see Fig. H lower panel). 

Turning to the InAs-based transistors, we note that 
different ratios of parameters describing SIA and BIA 
terms result in a different behavior of the conductance as 
a function of the transistor direction. For InAs the spin- 
orbit parameters are about order of magnitude larger, 
however, the electron mobility in InAs is low comparing 
to GaAs so the length of the quantum wire have to be 
shorter, otherwise we cannot assume ballistic transport 
regime. In calculations we have assumed L=4.0 /Ltmiii 
The results for the optimal angle (^=0° are presented in 
Fig- El They are compared to the case <^=135° which is 
the optimal angle determined within the linear approach. 
In Fig. we show conductance for opposite (upper panel) 
and parallel (lower panel) source-drain polarizations, re- 
spectively. 

The angle dependence of conductance for GaAs and 
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FIG. 4: Upper panel - Transconductance as a function of 
channel's direction. Lower panel - Relative difference of con- 
ductances for opposite and parallel source-drain polarizations. 



InAs devices is summarized in the upper panel of Fig. 0) 
where we plot the transconductance, IdG/dEpl, a quan- 
tity of particular experimental interest. The transcon- 
ductance in Fig. 01 is defined as the absolute value the 
derivative of conductance with respect to the Fermi en- 
ergy. For a given angle (p, we take into account the en- 
ergy region, where the conductance decreases. Indeed, a 
decrease of G is directly connected with the Datta-Das 
effect while its increases is connected also with the in- 
creasing number of modes participating in the transport. 
In the lower panel of Fig. 01 we present the Fermi energy 



dependence of another quantity which is also important 
in the experimental search for the Datta-Das effect,^ 
P = (G+_ - G++)/(G+_ + G++), where G++ and G+_ 
are conductances for parallel and opposite source-drain 
polarizations respectively. Our results imply, therefore, 
that there is much room for the enhancement of the 
Datta-Das signal by selecting an appropriate crystallo- 
graphic orientation of the device. 

In conclusion, we have developed a theoretical ap- 
proach suitable to numerical studies of baUistic elec- 
tron transport in semiconductor wires taking into ac- 
count spin-orbit effects. The model has been applied to 
the Datta-Das spin-FET of zinc-blende semiconductors. 
There are two main results of our model calculations. 
First, there is a strong dependence of the transistor con- 
ductance on its direction with respect to crystallographic 
axes. Secondly, our calculations explicitly show the im- 
portance of nonlinear terms in BIA Hamiltonian. For 
particular model of the spin-orbit parameters, Eq. (01 
and Eq. (j^J , we have also determined transistor direction 
optimal for the observation of the Datta-Das effect, which 
is [120] and [100] for the GaAs and InSb-based devices, 
respectively. 
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